Quantum Quench of an Atomic Mott Insulator 
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We study quenches across the Bose-Hubbard Mott-insulator-to-superfluid quantum phase transi- 
tion using an ultra-cold atomic gas trapped in an optical lattice. Quenching from the Mott insulator 
to superfluid phase is accomplished by continuously tuning the ratio of Hubbard tunneling to interac- 
tion energy. Excitations of the condensate formed after the quench are measured using time-of-flight 
imaging. We observe that the degree of excitation is proportional to the fraction of atoms that cross 
the phase boundary, and that the quantity of excitations and energy produced during the quench 
have a power-law dependence on the quench rate. These phenomena suggest an excitation process 
analogous to the Kibble-Zurek (KZ) mechanism for defect generation in non-equilibrium classical 
phase transitions. 
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The study of non-equihbrium phase transitions is wide 
ranging, touching on topics as diverse as the formation 
of structures in the early universe [1] and the practical- 
ity of adiabatic quantum computing [2]. The so-called 
"Kibble-Zurek" (KZ) mechanism has been used to un- 
derstand some universal features — principally the rate of 
topological defect formation — of quenches across classical 
phase transitions [l], Q . "Quench" in this context refers 
to rapidly varying a thermodynamic parameter in order 
to drive the system out of equilibrium. The KZ theory 
has recently been extended to quantum phase transitions 
[4-15|. In contrast to the classical case, quantum phase 
transitions involve completely closed quantum mechan- 
ical evolution at zero temperature, for which quenches 
are accomplished by varying a parameter in the Hamilto- 
nian in order to tune between different quantum phases. 
While the KZ mechanism has successfully been tested 
for classical transitions (e.g., on liquid crystals [3]), and 
spontaneous vortex formation has been observed during 
cooling an atomic gas through the Bose-Einstein conden- 
sation transition [17] , experimental examination of quan- 
tum quenches have been scant. Notably, there is evidence 
that the formation of ferromagnetic domains in a spin-1 
Bose-Einstein condensate can be attributed to a sudden 
quantum quench In this work, we probe quantum 

quenches for a paradigm of quantum phase transitions — 
the Bose-Hubbard (BH) model — using atoms confined in 
an optical lattice. In contrast to previous experiments 



20l|, we quench from the Mott-insulator (MI) to the 
superfluid (SF) state, and we systematically investigate 
the formation of excitations as the quench amplitude and 
rate are varied. 

In our experiment, a cubic optical lattice formed from 
three intersecting pairs of 812 nm laser beams is super- 
imposed on a parabolically confined ^^Rb Bose-Einstein 
condensate; details of our apparatus can be found in Ref. 
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and references therein. The atoms in the lattice are 
described by the inhomogeneous BH model with tunnel- 
ing energy t and interaction energy U , the ratio of which 



is controlled by tuning the lattice laser intensity to ad- 
just the lattice potential depth s. By changing s, MI and 
SF phases can be sampled inhomogeneously in the gas 
[2^ : for s > ISEji, nested Mott-insulator and superfluid 
layers exist in the lattice, and for s < ISiJ/j the gas is 
purely superfluid, as shown in Fig. 1 {En = h/8m(fi, 
where m is the atomic mass, d — 406 nm is the lattice 
spacing, and h is Planck's constant). 

Quenching across the SF-MI phase transition is ac- 
complished by adjusting s dynamically in such a way 
to transform the gas between equilibrium configurations 
with and without atoms in the MI phase present (Fig. 1). 
While quenches are possible on all relevant timescales, in 
this paper we explore quenches that occur at rates 1/tq 
that are too slow to excite atoms into higher vibrational 
states in the lattice potential. How 1/tq compares with 
the Hubbard energies is complicated because the phase 
boundary is crossed at a range of densities and t/U in 
the trap, and therefore the Hubbard energies Uc and tc at 
the phase transition change during the quench. Despite 
this, the quenching rate is always slow compared with 
Uc- 1/tq varies from 1 x 10^''-0.2 Uc/h. The quench 
rate is not consistently fast or slow compared with tc [l^ 
or the confining trap frequencies, the geometric mean of 
which varies from 43 ± 2 Hz at s = OEn to 82 ± 6 Hz at 
s = 25Er. 

The effect of variations in the fraction of atoms crossing 
the SF-MI phase boundary is investigated by quenching 
s linearly in 5 ms to s = AEn (corresponding to t/U ~ 1, 
i.e.,_the purely SF regime), as shown in the inset to Fig 
2 24;]. The ratio t/U changes non-linearly during this 
quench; in the large s limit, t/U oc e~^^. The fraction 
of atoms crossing the phase boundary is varied by adjust- 
ing the initial lattice potential depth sq- The data shown 
in Fig. 2 sample a range such that at high sq nearly all 
of the atoms start in the MI phase (with fillings ranging 
from 1-3 particles per site) and therefore cross the phase 
boundary; and at low enough sq so that all of the atoms 
are in the SF phase, and consequently no atoms cross the 
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FIG. 1. Quench across MI-SF boundary (top) and excita- 
tion measure (bottom). As shown by the vertical lines, the 
trapped gas samples a range of densities and effective chemi- 
cal potentials p. (in the LDA) [S^]; MI regions in the gas are 
colored red and SF blue. Given the overall confining potential 
in our experiment and atom number ((161 ± 13) x 10'^, av- 
eraged across all measurements), the maximum ft is roughly 
fixed at 2U for the measurements described here, which cor- 
responds to a central filling of 4 atoms per site at low s and 
a central MI with 3 atoms per site at high s. A quench is ac- 
complished by rapidly reducing s and thereby increasing t/U, 
as shown by the arrow. A slice (black line) through a typical 
image (inset) taken after a quench for so = 25Er is displayed. 
The image is fit to a smooth profile (red line), which is used 
to determine the deviation Xij (blue line) at each pixel in a 
masked region (gray). 



phase boundary. The fraction of atoms in the MI phase 
before the quench, which is identical to the overah frac- 
tion of atoms traversing the phase boundary, is shown 
as a red Hne in Fig. 2, and is determined according 
to the measured atom number and a zero-temperature 
mean- field calculation in the LDA [g^l . After the quench, 
the lattice is turned off in 200 /iS. This "bandmapping" 
step [2l'| — which maps quasimomentum in the lattice to 



free momentum and suppresses atom diffraction — is nec- 
essary to improve the imaging signal-to-noise ratio given 
the long expansion times employed for these measure- 
ments. 

The amount of excitation produced during the quench 
is determined by measuring the deviation from a smooth 
profile of time-of-flight images taken after release of the 
trapped gas immediately following bandmapping. The 
absorption image is taken after a relatively long 50 ms 
of f ree expansion, so that vortices, if present, are visi- 
ble 
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and phase gradients related to other topological 
or wave-like excitations are converted into large density 

H-IHi. 



fluctuations I26H281. We fit the image to a smooth func- 
tion / that is a combination of a Thomas-Fermi pro- 
file and a Gaussian, and measure the amount of ex- 
citation as the deviation from the smooth profile: 



lY^ij ODij, where i,j 



■It] " fij I 

index the pixels in the image within a mask defined by an 

imaging signal-to-noise ratio greater than 5, OD is the 
measured optical depth, and a is a proportionality con- 
stant that is determined using a numerical simulation. 
We find that all of the images used in this work are well 
described by this fit — a condensate appears present after 
the quench under all circumstances, and the condensate 
fraction varies from 0.35-0.6 across all of the data. While 
it was suggested in Ref. [2^ that the condensate fraction 
may oscillate after the quench, we find no evidence for 
such behavior. 

The measure is chosen such that it is related to 
the fraction of atoms in excited states for the trapped, 
weakly interacting gas present before bandmapping. The 
physical meaning of can be understood most straight- 
forwardly for a one-dimensional non-interacting gas. 
In this case, the density profile after sufficiently long 
TOF is the momentum distribution n{q) = \il>{q)\^ = 
l'0o('?) + <^'0('?)|^i where (^t/i are plane-wave excitations, 
i^ail) = -\/no((7) is the ground-state condensate wave- 
function, and we work in the momentum representa- 
tion. After averaging over random excitation phases, 
the number of atoms in excited states is ^ dq \5ij}(ci)\ — 
J dq [{n{q) ~ no{q)]'^ /2no{q). Given that J dq no{q) is 
the total number of atoms, x^ is naturally interpreted 
as proportional to the number of excited atoms in the 
non- interacting limit. 

Using a numerical simulation of the three-dimensional 
time-dependent Gross-Pitaevskii equation, we deter- 
mined both that x^ accurately reproduces the fraction of 
Bogoliubov excitations for a trapped condensate and the 
constant a. We start the simulation with a condensate at 
equilibrium in a parabolic trap (using the experimental 
atom number and trap parameters) , and imprint Bogoli- 
ubov excitations under the LDA for a range of wavevec- 
tors corresponding to 0.8 — 3 jim^^; the Thomas- Fermi 
radius of the gas is approximately 10 /im before release. 
Images are generated by time evolving the condensate 
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wavefunction for a free expansion and then integrating 
through the imaging hne-of-sight. The measure is 
determined for a range of excitation fractions averaged 
over 10 relative phases. We determine that is equal 
to the fraction of excited atoms for a = 10 under simu- 
lated conditions. This method is an approximation, and 
does not properly account for long-wavelength (i.e., trap- 
length-scale) excitations or topological excitations such 
as vortices, which are evident in the insets to Figs. 1 and 
2. 

As shown in Fig. 2, we find that the amount of excita- 
tion is proportional to the fraction of atoms crossing the 
phase boundary. Below the emergence of the unit filling 
MI phase at so ~ ISEr, is constant at Xo ~ 0-06 (de- 
termined by averaging over all images with sq < ISEji, 
and indicated by the dashed line in Fig. 2), a value that 
is consistent with the combination of photodetection shot 
noise and technical noise present in our imaging system. 
Above So ~ 13i?j^, the degree of excitation grows, until 

saturates to approximately 0.17 at high enough lattice 
depth, for which more than 90% of the atoms are in the 
MI phase. 

The behavior evident in Fig. 2 suggests that a Kibble- 
Zurck-like mechanism is responsible for generating exci- 
tations during the quench. In the KZ picture, the diverg- 
ing relaxation time near the phase boundary "freezes in" 
fluctuations in the relative phase between atomic wave- 
functions at different lattice sites present in the MI Q. 
Some time after crossing the phase boundary, dynam- 
ics effectively restart, and the fluctuations develop into 
superfluid excitations, potentially including sound waves 
and topological excitations such as vortices. Given that 
only the regions of the lattice that cross the SF-MI phase 
boundary will give rise to excitations, the direct relation 
between the fraction initially in the MI phase and the 
degree of excitation is strong evidence for KZ physics. 

In the KZ scenario, the quench rate controls the num- 
ber of excitations generated according to a power law 
that depends on the critical exponents for the phase tran- 
sition. We measured this power law, as shown in Fig. 3, 
across two orders of magnitude in quench rate. For this 
measurement we quench the lattice potential depth start- 
ing from a gas composed nearly entirely of the MI phase 
at So = 20Er (i.e., t/U = 0.005) according to s(r) = 

0.25 In^ + e~2v^) (a^ w 5 nm is the scattering 

length) so that the rate of change d(t/U)/dT w I/tq is 
approximately constant [s^. In addition to measuring 
how depends on we also measure the kinetic 

energy generated by the quench. The kinetic energy per 
particle KE is measured from TOF images according to 
KE = TO (r^) /2T(gj, with the second moment of the 
density distribution (r^) = 3/2 • Y.^J OD,Jrl/J2^J OD,j 
after the TOF ttof- The factor of 3/2 arises from as- 
suming the energy is distributed equally among three di- 
rections. 
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FIG. 2. Amount of excitation produced by quenching as the 
fraction crossing the phase SF-MI boundary is varied. The 
top left inset shows the experimental timeline; a magnetic 
field gradient is applied to support the atoms against gravity 
during TOF. Characteristic TOF images are shown as insets 
for So =9, 15, and 25 Er, and the error bars in this and 
the next figure are the standard deviation for the average 
taken over 5 images. There is a 7% systematic uncertainty to 
So, which is calibrated using Kapitza-Dirac diffraction. The 
overall uncertainty in the MI fraction (red line) ranges from 
30% at So = 16Er to 10% at so = 20Er- below so = 12Er 
and above so = 22Er the uncertainty in MI fraction is zero. 



A fit to the data in Fig. 3 reveals power laws 1/tq 
for x^ and KE consistent within the fit uncertainty: r — 
0.31 ± 0.03 and 0.32 ± 0.02, respectively. While there are 
numerous detailed theoretical predictions for the number 
of excitations produced during a quench across the MI- 
SF phase transition [1, S [H - 14 1 , none that we know of 
are directly applicable to our experiment. Generically, 
the size of domains associated with an excitation formed 
during the quench should scale as j-q^'^^'^'^'^ [31|, where 
z is the dynamical critical exponent and the correlation 
length diverges as ^ ^ {t/U — tc/Uc)'" near the phase 
boundary. In three dimensions, the density of excita- 
tions is therefore propor tional to I/t^^'^^'+^I For our 
experiment, nearly all of the atoms cross the "generic" 
phase transition and not the multi-critical point at the 
the "tip" of the MI "lobes." In this case v ~ \/2 and 
z = 2, and therefore the number of excitations should 
scale as I/tq'^, which is inconsistent with our data. This 
disagreement may be explained by numerous issues that 
deserve more theoretical attention. For example, the 
spatially inhomogeneous nature of the phase transition 
that gives rise to a phase transition "front" that moves 
through the gas has been examined in the context of the 
classical SF phase transition [s^l, but not for the quan- 
tum case. Also, the finite size of the gas will affect quench 
dynamics, as discussed in Ref. 15| for the BH model in 
ID. Finally, since the data here were taken at low but 
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FIG. 3. Excitation dependence on quench rate. The amount 
of excitation and kinetic energy determined from TOF images 
is shown as the quench rate 1/tq is varied. The measured 
offset Xo is subtracted from x^. Analogously, the measured 
expansion energy without the quench KEo (determined by 
averaging across images with so < 13Er) is subtracted from 
KE. The inset shows the experimental timeline, and the solid 
lines are power-law fits to the data. 

finite temperature (the initial condensate fraction was 
more than 90% before turning on the lattice), thermal 
effects may play an important role in the quench dy- 
namics; see Ref. for a analysis in the context of the 
Sine-Gordon model. 

In conclusion, the method we have demonstrated pro- 
vides a window into excited states and dynamics, which 
are beyond our current theoretical understanding in a 
wide variety of strongly interacting many-body quan- 
tum systems. Quench dynamics may also have signifi- 
cant consequences for thermometry in optical lattice ex- 
periments [33|. One commonly employed technique to 
estimate temperature in a lattice is to slowly turn off 
the lattice potential, measure temperature, and then in- 
fer entropy in the lattice assuming that the turn off was 
adiabatic. We find across a wide range of linear lattice 
quench rates that adiabaticity is violated; for example, 
for a quench from sq = 20£'i{, decreases from 0.17 
to only 0.12 for turn off times varying from 5 to 25 ms. 
Finally, in the future, optical lattice quenches may be in- 
vestigated in other contexts, such as reduced dimensions 



and fermionic gases. 

The authors acknowledge financial support from 
the National Science Foundation (award 0448354) and 
DARPA OLE program. Any opinions, findings, and con- 
clusions or recommendations expressed in this material 
are those of the authors and do not necessarily reflect the 
views of the National Science Foundation. 



* Now at: Stanford Research Systems Inc., 1290-D Ream- 
wood Ave, Sunnyvale, CA 94089 
^ Now at: Department of Atmospheric Sciences, University 
of Alaska, 930 Koyukuk Dr, Fairbanks, AK 99775 

[1] T. Kibble, J. Phys. A 9, 1387 ( 1976). 

[2] R. Schiitzhold et ai, Phys. Rev. A 74, 060304 (2006) 
[3] W. Zurek, Nature 317, 505 (1985). 
[4] F. Cucchietti et ai, Phys. Rev. A 75, 23603 (2007) 
[5] J. Dziarmaga, Adv. Phys. 59, 1063 (201 0). 
[6] J. Dziarmaga et al, Phys. Rev. Lett. 101, 115701 (2008)| 
[7] C. De Grandi et al, Phys. Rev. B 81, 224301 (2010)J 
[8] C. DeGrandi and A. Polkovnikov, Quantum Quenching, 
Annealing and Computation, edited by A. K. Chandra 
et al.. Lecture Notes in Physics 802, 75-114 (Springer 
2010). 

[9] L. Mathey et al, Phys. Rev. A 81, 033605 (2010)| 
[10] L. Mathey et al, Phys. Rev. A 80, 041601 (2009)| 
[11] M. Schiro et al, ar Xiv:1102.1658 (2011). 
[12] R. Schiitzhold et al, Phys. Rev. Lett . 97, 200601 (2006) 
[13] M. Snoek, arXiv:1101.5006 (2011). 
[14] C. Trefzger et al, arXiv: 1008. 1285 (2010). 
[15] G. Roux, Phys. Rev. A 81, 053604 (2010) 
[16] I. Chuang et al. Science 251, 1336 (199ljj 
[17] C. Weiler et al. Nature 455, 948 (2008) 
[18] L. Sadler et al. Nature 443, 312 (2006) ,^ 
[19] M. Greiner et al. Nature 419, 51 (2002)] 
[20] W. Bakr et al. Science 329, 547 (2010) 
[21] D. McKay et al, Phys. Rev. A 79, 063605 (2009)' 
[22] B. DeMarco et al, Phys. Rev. A 71, 063601 (2005)1 
[23] For our fastest quenches, 1/tq is always fast compared 
with tc/h (1/tq: 10-100 tc/h); for the slowest quenches, 
l/rg is comparable to or slower than tc/h {1/tq: 0.01- 
1 tc/h); and for intermediate speeds, 1/tq may be fast 
compared with tc/h at the beginning of the quench (e.g., 
at s = 25-E_r) and slow at the end. 
[24] A 5 ms total quench time was chosen for all data in this 
work so that the finite excitation (exponential) lifetime — 
measured to be ~ 30 ms after a quench from so — 2QEb. — 
had little effect. 

[25] K. Madison et al, Phys. Rev. Lett. 84, 806 (2000)j 

[26] S. Dettmer et al, Phy s. Rev. L ett. 87, 160406 (200T)1, 
[27] D. Petrov et al, Phys . Rev. Le tt. 87, 050404 (2001) 
[28] D. Hellweg et al, Appl. Phys. B-Lasers O. 73, 781 (2001). 
[29] E. Ahman et al, Phys. Rev. Lett. 89, 250404 (2002) 
[30] W. Zwerger, J. Opt. B: Quantum Semiclass. Opt 5, S9 
(2003). 

[31] A. Polkovnikov, Phys. Rev. B 72, 161201 (2005). 
[32] A. del Campo et al, arXiv:1010.6190 (2010). 
[33] D. McKay et al, arXiv: 1010.0198 (2010). 



